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A study of large field configurations in MC simulations 
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We discuss a new approach of scalar field theory where the small field contributions are treated perturbatively 
and the large field configurations (which are responsible for the asymptotic behavior of the perturbative series) are 
neglected. In two Borel summable A<^* problems improved perturbative series can be obtained by this procedure. 
The modified series converge towards values exponentially close to the exact ones. For A larger than some critical 
value, the method outperforms Fade's approximants and Borel summations. The method can also be used for 
series which are not Borel summable such as the double-well potential series and provide a perturbative approach 
of the instanton contribution. Semi-classical methods can be used to calculate the modified Feynman rules, 
estimate the error and optimize the field cutoff. We discuss Monte Carlo simulations in one and two dimensions 
which support the hypothesis of dilution of large field configurations used in these semi-classical calculations. We 
show that Monte Carlo methods can be used to calculate the modified perturbative series. 



In a recent publication Q , we have shown with 
three nontrivial 0"* problems that cutting off the 
integration range of </) at each site modifies the 
large order behavior of these series. The modi- 
fied series converge toward values which are ex- 
ponentially close to the exact ones. The ex- 
amples were the anharmonic oscillator and the 
Landau-Ginzburg hierarchical model which are 
Borel summable, but also the double- well which is 
not Borel summable, and in this case, the method 
incorporates instanton effects. More generally, we 
are interested in developing approximate treat- 
ments where the small field configurations are 
treated perturbatively and the large field config- 
urations semi-classically. 

All the examples discussed above can be solved 
accurately with numerical methods H,^. In the 
following, we discuss the possibility of using the 
Monte Carlo method to perform such calcula- 
tions. We consider scalar models with one com- 
ponents at each site. We define the norm of 
a configuration C as |C| — MaXxW<i>x^- The 
norm distributions of MC configurations for A(/)^ 
models in 1 and 2 dimensions are shown in Fig. 



1. One sees that the distributions have rapidly 
falling tails and that the number of large norm 
configurations is exponentially suppressed. One 
would thus expect that the error on correlation 
functions caused by cutting-off the configurations 
with |C| > (^raax has the generic form 



Error - 



(1) 



In addition, the sites with large fields are clus- 
tered, and the configuration has a typical shape 
in these regions as shown in Fig. 2. This suggests 
that the large field configurations contribute di- 
lutely to the functional integral. This is just what 
is needed for a semi-classical calculation. 

The similarity of the (Borel summable) results 
found in Ref. suggests that in general the cor- 
rections due to the field cutoffs can be expressed 
as simple one-dimensional integrals. The correc- 
tion to the zeroth order of the ground state of the 
anharmonic oscillator (in other words, the correc- 
tion to the ground state of the harmonic oscilla- 
tor) has the semi-empirical form [Q. 
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Figure 1. Norm distribution for X<p'^ models in 1 
and 2 dimensions 



Figure 2. A large norm configuration for the an- 
harmonic oscillator 
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Fig. 3 shows that this is a good approximation 
which becomes better as 4>max increases. 

When (pmax is large, the semi-empirical formula 
has the asymptotic form 

SE^°^ ~ 2^-i/20™,,e-*"- (3) 

This correction has the same functional depen- 
dence as the semi-classical splitting for the first 
two levels of the double- well jj] : 

^4°)oc(|)^/^e-" (4) 

provided that we take 5*0 — (t^inax^ the clas- 
sical action from the classical configuration 
'/'maxe~'^~^°' . This confirms the validity of the 
dilute gas approximation. Calculations of higher 
order perturbative coefficients and in higher di- 
mensions are in progress. 



Figure 3. Ratio of the estimate and the numerical 
value for Se'^^^ vs (t)max- 

It is possible to use the MC method to calcu- 
late the coefficients of the modified perturbative 
series. This will be particularly useful to check 
the validity of semi-classical methods in higher 
dimensions. In the case of the anharmonic os- 
cillator H = + 9^*^/2 + A^'', the first order 
coefficient is 

< ^ (/.^ >Gaussian / ^ = 3/4 . (5) 

X X 

We have calculated the modified averages corre- 
sponding to various field cuts. The results are 
shown in Fig. 4 and are in good agreement with 
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the accurate numerical results. Note that these 
coefficients have a sizable lattice spacing depen- 
dence which can taken care with various methods 
which will be discussed elsewhere 
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Figure 4. First perturbative coefficient for the 
ground state of the anharmonic oscillator with 
different field cuts : 1 is 1.5, 2 is 1.75, 3 is 2 etc. 
. The continuous line is the numerical answer. 



At a fixed value of (fimax, the error due to the 
field cut decreases when the coupling increases. 
At a fixed order in the modified perturbative se- 
ries corresponding to the same fixed (j)max, the 
error due to the omission of the higher orders 
decreases when the coupling decreases. Conse- 
quently, at fixed ipmax, there is a value of the 
coupling where the the accuracy is optimal. This 
value is approximately where the two error curves 
meet. Conversely, at fixed coupling, one can find 
an optimal value of (pmax ■ This value can be re- 
lated to the optimal values of variational param- 
eters used in Refs. In these calculations, 
an intermediate mass term is introduced. The ef- 
fect of this term is at the same time to reduce 
the contributions of the large fields and to re- 
duce the argument of the exponential which is 
expanded. Simple power law relations appear 
in the strong and weak coupling limits for the 
problem of the one dimensional integral ^ . 

The coupling dependence of the optimal (l)max 
is a question which needs to be understood if we 
want to apply similar methods to gauge theories. 



In Wilson's formulation, at fixed lattice spacing, 
a field cut inversely proportional to the gauge 
coupling is imposed when one uses the compact 
group integration. Consequently, the possibility 
for scalar models of having an optimal (jimax vary- 
ing like the inverse square root of the quartic cou- 
pling in some regime would be quite interesting. 
We have started investigating the possibility of in- 
troducing field cuts in gauge theories by using the 
pubhc OSU quenched configurations [||. These 
configurations are in the Landau gauge, and the 
distributions of values of a fixed matrix element 
have similarities with the scalar case. 

Finally, we would like to mention that, in the 
large- iV limit, the critical potentials of the 0{N) 
models in 3-dimensions have a finite radius of con- 
vergence when expanded in the 0{N) invariant 
quantity 0^ This is due to two conjugated 

branch cuts in the complex 0^ plane. It is tempt- 
ing to consider the model with a field cut which 
corresponds to this radius of convergence. How- 
ever, this procedure generates sizable errors com- 
pared to the reliable procedure which consists in 
defining the critical potential using Fade approx- 
imants. 
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